Magnetic multipole analysis of kagome and artificial ice dipolar arrays 
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We analyse an array of linearly extended monodomain dipoles forming square and kagome lattices. 
We find that its phase diagram contains two (distinct) finite-entropy kagome ice regimes — one disor- 
dered, one algebraic — as well as a low-temperature ordered phase. In the limit of the islands almost 
touching, we find a staircase of corresponding entropy plateaux, which is analytically captured by 
a theory based on magnetic charges. For the case of a modified square ice array, we show that the 
charges ('monopoles') are excitations experiencing two distinct Coulomb interactions: a magnetic 
'three-dimensional' one as well as a logarithmic 'two dimensional' one of entropic origin. 
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An intensifying experimental effort is under way 
aimed at constructing dipolar arrays of monodomain 
nanomagnetsii"— From a technological viewpoint, this is a 
promising approach as it harnesses mature techniques — 
for example those involved in the lithographic prepara- 
tion of arrays — often highly developed in an industrial 
context. Indeed, one of the fundamental questions from 
this perspective is to understand interactions and dynam- 
ics of closely packed nanomagnets as these are among the 
factors determining performance limits on high-density 
magnetic storage media. 

From a many-body physics angle, this forms one as- 
pect of a wider program of studying artificial systems 
which can be constructed, manipulated and probed on 
the level of an individual degree of freedom. Pioneer- 
ing experiments by Wang et al. on the two-dimensional 
'square ice' array have established a proof-of-principle of 
the feasibility of such an enterprise^ by demonstrating 
that the magnets assume configurations reflecting their 
interactions, paving the way to a study of their collective 
propertiesiii^ 

While the study of dipolar magnets is a rich subject 
with a venerable tradition^ with a recently rejuvenated 
activity in the context of cold gasea^ — interesting on ac- 
count of the long ranged and anisotropic nature of the 
interactions — recent work on the spin ice compounds^ 
has found a number of surprising new phenomena. In 
particular, it was found that even the complex dipolar 
interactions can lead to large low-temperature degenera- 
cies, apparently with little additional fine-tuning.—"— A 
simple model which considered each magnetic dipole as a 
pair of nearly opposite magnetic charges ( "monopoles" ) 
was able to account for this phenomenon analytically, 
and established that in spin ice, these monopoles are the 
appropriate degrees of freedom at low temperatures 

In this paper, we analyse interactions in artificial dipo- 
lar arrays in d = 2. We study two geometries, the qual- 
itative behavior of which varies reflecting simple lattice 
parameters such as even/oddness of their coordination. 
Besides the abovementioned artificial spin ice, we partic- 
ularly analyse the case where the dipolar islands reside 
on the bonds of a honeycomb lattice with their centres 
forming a kagome array, which has been a focus of re- 
cent experimental studiesi^"— For this kagome array. 



we find a series of regimes: firstly, the paramagnct cor- 
responds to a gas of charges ±1 and ±3. As the inter- 
action is increased in strength, there is a crossover to a 
gas of charges ±1 (kagome ice I, Kl) which terminates in 
an ordering transition of the charges in a two-sublattice 
NaCl structure on the honeycomb lattice (kagome ice II, 
K2). There being exponentially many configurations of 
the islands which yield a NaCl charge configuration a fi- 
nal transition into an ordered structure with tripled unit 
cell removes this entropy. In a particular limit, we can 
extend the monopole theory for spin ice to a multipole 
one which captures all these regimes in analytic detail. 

We supplement this by an analysis of the square ice ar- 
ray, where we find a hybrid Coulomb interaction between 
monopoles: one is the usual 1/r interaction between mag- 
netic charges, whereas the other is of entropic origin and 
exhibits a logarithmic distance dependence. An addi- 
tional linear confining force can be removed by shifting 
the monopoles in the third dimension so that their loca- 
tions define tetrahedra. In closing, we remark on connec- 
tions to other systems as well as issues of dynamics and 
finite-size effects. 

Dipolar needles: The basic degree of freedom is a 
dipolar island which we model with the following (ide- 
alised) properties. The magnetic moment density pi/l 
is uniform along the length, /, of the island. The is- 
lands are monodomain with the moment pointing along 
their axis, yielding an effective Ising spin — as was verified 
experimentallyii"— ii^rJ^ We take the shape of the island 
to be that of a needle, with a length much greater than 
its width, neglecting its finite transverse extent. A nee- 
dle with uniform magnetisation density q= sets up 
the same potential U{R) as that of a pair of magnetic 
charges ±q located at its tips (in the context of micro- 
magnetic domain walls, see Refs. The value of 
the Ising spin encodes which of the ends — located at 
and ff, — hosts the positive charge: 
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In the following, we express lengths in units of a, the 
nearest-neighbor distance of the respective lattices. The 
dipolar islands are located on the bonds of the lattice. 



2 



In the regime l/a ^ 1, which is approximately rcahzcd 
in recent experiments?^^— it is convenient to define e = 
1 — I /a. We generally use q^/a as the unit of energy, 
but shall also refer to the nearest neighbor Ising coupling 
strength of the islands, J, where this is more instructive. 

Many of the salient features of these spin systems can 
be deduced intuitively in the charge picture which, in the 
limit of e — >■ 0, admits an analytic treatment. 

Let us first consider the case of the kagome array. 
The vertices of the kagome lattice can have total charge 
Q = ±1 (one-out, two-in or two-out, one-in) or Q = ±3 
(all in or all out). The corresponding on-site energies of 
El = and £^3 = are such that the latter appear 
only at high temperature, T, in a phase that can be clas- 
sified as a paramagnet. Below TJ^^ ~ £'3 — i?i ~ 2 J, the 
population of charge ±3 vertices is exponentially sup- 
pressed. This crossover occurs without change of sym- 
metry. The remaining configurations with Q = ±1 ev- 
erywhere are precisely those satisfying the ice-rule (two- 
in, one-out or vice versa). The tendency to this kind 
of "order" was investigated and confirmed in recent 
experiments i^— 

Unlike for the case of spin-ice, where a coordination 
number of 4 allows for vertex charges to vanish (and fi- 
nite charges arc excitations in d = 3, see Ref. Il2i ). the 
odd coordination of 3 enforces the presence of magnetic 
charges in kagome ice at any temperature. 

This kagome ice (Kl) regime is a highly frustrated 
phase of matter, with an entropy (known exactly) about 
5/7 that of a free spin, S^^ « § In « 0.501. 

At lower temperatures T;^^ ~ (f'/o-, interactions be- 
tween vertices assert themselves. They lead to a "charge- 
ordered" NaCl ionic crystal, in which each sublattice of 
the honeycomb lattice only hosts charges \Q\ — 1 of a 
given sign — the resulting magnetic charge order is thus 
described by an Ising order parameter. 

While the NaCl structure is fully ordered in terms of 
charges, it is still consistent with an exponentially large 
number of configurations of the underlying Ising degrees 
of freedom of the islands. This kagome ice II (K2) regime 
is therefore only partially ordered. To calculate this en- 
tropy we map each spin configuration in kagome ice II 
to a dimer covering of the hexagonal lattice: note that 
at any vertex, there is precisely one minority magnetic 
charge. As Qa = —Qp, the edge with the minority charge 
is necessarily connected to the minority charge at the 
neighboring vertex, as shown in Fig. [1^). Thus, each 
charge configuration singles out an (oriented) dimer cov- 
ering, with precisely one dimer pointing from the neg- 
ative to the positive charge impinging on each vertex. 
The orientation of the dimcrs encodes the broken Ising 
symmetry. 

The entropy of hexagonal dimer coverings is known by 
a mapping to the triangular Ising antiferromagnet with a 
well-known groundstatc entropy ^^pj^. Adjusting for the 
number of elementary degrees of freedom — islands rather 




FIG. 1: (color online) a) Orientations of Ising spins on the 
links are mapped to a distribution of magnetic charges lo- 
cated at the extremities of each island. In kagome-ice, the 
total charge at each vertex is constrained to \Q\ = 1. At low 
temperatures, vertex charges order such that all vertices of 
sublattice a have charge Qa = (—1)" (as shown). To count 
the configurations corresponding to that NaCl charge-state, 
we map these to dimer coverings of the honeycomb lattice, 
where oriented dimers are present at all links connecting the 
minority charges of the adjacent vertices, b) To count the 
dimer coverings of the hexagonal lattice, we note that dimers 
present on the hexagonal lattice can be mapped to broken 
links of the antiferromagnet on the dual triangular lattice. 
The particular configuration shown contains the unit-cell of 
the groundstate, formed by the three highlighted plaquettes. 
Arrows shown at the vertices display the direction of the local 
dipole moment. 



then dimers — we find an entropy per island of 

5iii = ^5^FM« 0.108. (1) 

Finally, at very low temperatures, further neighbor 
interactions induce a transition into an ordered phase 
which removes all remaining entropy. As the most 
favourable arrangement of a single hexagonal plaquette is 
for all islands to point head-to-tail,— it is natural to guess 
that the state should optimise the occurrence of that mo- 
tif. However, no full tiling of the honeycomb lattice with 
such oriented loops is possible. Instead, one finds a con- 
figuration, in which 2/3 of all plaquettes are in such a 
state, which is precisely the one shown in Fig. [T] The 
three plaquettes highlighted in Fig. [T}d) form the tripled 
X \/3 unit-cell of the crystal with periodicity in the 
lattice vectors vi — (0,3a) and V2 = (3-\/3a/2, 3a/2). We 
will refer to this groundstate as the 'loop state', below. 
Note that the plaquettes without an oriented loop host 
no dimers at all. 

Analytic magnetic multiple expansion: The cascade 
of regimes outlined here can analytically be accounted 
for in the limit e — > 0, as each transition/crossover is 
associated with a parametrically different energy scale. 
This approach amounts to a development of a multi- 
pole expansion beyond the dumbbell model of spin ice. 
One finds for the energy of the loop state (expressed per 
hexagonal site): 

^ g2 r _2 3 a, 3 , S3, .I 
^-^{7!^ + (2-2) + r+(^+4 + 2)^}- 
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The 0{l/e) term enforces the |Q| = ±1 constraint; it is 
due to the interactions between the charges on the same 
vertex. At next order, the honeycomb Madelung energy, 
with a « 1.5422197 drives the Kl K2 Ising transition. 
This term reflects the long-range 1/r interaction between 
the charges (monopoles) of different vertices. The 0{e^) 
term docs not distinguish between dimer configurations 
and hence docs not lift their degeneracy; this instead 
happens at O(e^), where the dipole-dipole term, with 7 r:! 
—2.226947, rather than the monopole-quadrupole term 
with 5 « —0.5829489, leads to a lifting of the degeneracy 
between dimer states. (In passing, we note that these 
constants can be accurately evaluated using finite-size 
scaling converging as 1/L'^ by arranging the unit cell of 
the summation to have vanishing monopolc, dipolc and 
quadrupole moments). 

The order of this term is easily understood. The 
charges ±q are displaced a distance ^ ea from the hexag- 
onal vertex they belong to, which generates a dipole mo- 
ment of size \p\ ~ eaq and hence an interaction energy of 
0{e^). This dipole moment points in one of three direc- 
tions, namely along the dimer emanating from the site, 
so that sites connected by a dimer have aligned dipole 
moments. The loop pattern optimises the local dipo- 
lar energy, in that each dimer is oriented favourably with 
respect to its four nearest neighbours and due to the con- 
vergent nature of the 1/r^ dipolar interactions in d = 2, 
there is little scope for this state being destabilised. 
Numerics: Similarly, in numerics, the regimes are well 
resolved for small e, where they correspond to extended 
plateaux in thermodynamic quantities. The data shown 
here are for e = 0.05 from Monte Carlo simulations on a 
lattice of linear size L = 24, with the interaction summed 
over periodic copies up to distance A = 500 L. Wc em- 
ploy a mixture of single-spin flips and the short loop 
algorithm^i^ which seeks a chain of aligned spins and de- 
fines a loop when this chain first crosses itself. As the 
presence of different order parameters nonetheless ren- 
ders a fully ergodic algorithm difficult to achieve we have 
pinned down the final ordering transitions by heating and 
cooling runs, the latter at times leading to the formation 
of large-scale domains. 

We find a broad paramagnet — )- Kl crossover, followed 
by a continuous and hysteresis-free onset of the Ising 
charge order (Kl — > K2). These phases exhibit the pre- 
dicted entropy, which drops to zero at T = Torder, below 
which we find the loop crystal, the six-fold degeneracy 
of which is due to the (three-fold) translational symme- 
try breaking, on top of the two-fold Ising spin reversal 
symmetry. 

Interestingly, the orientation of the dipolcs for a given 
translation of the crystal is determined in the transition 
at the higher temperature T-^J^. As this ultimately deter- 
mines the orientation of the loops, one could also refer 
to it as a chiral transition. Further, the acceptance rate 
of the loop moves of our Monte-Carlo simulations ap- 
proaches unity in the kagome ice K2 regime, indicating 
that it can be regarded as a loop gas. 




e=l-l/a T [Q^/a] 



FIG. 2: (color online) Numerical results for the dipolar 
kagome lattice: a) phase-diagram locating the phases (from 
the peaks in the heat capacity) as a function of the length of 
dipolar magnetic islands e = 1 — l/a and temperature in units 
of the NN-interaction Ji. The inset details the behaviour 
for small e: quadratic scaling of the ordering temperature as 
Torder ~ 9^/" ^s predicted by Eq. ([2]), and a constant tran- 
sition temperature T^J. ~ 9^ /a- b) An exemplary trace of 
the heat-capacity C and entropy per spin S as a function of 
temperature for e = 0.05. The transitions are observed as 
well separated peaks in the heat capacity. The entropy shows 
broad plateaux matching the theoretical predictions (dotted 
lines) for the different phases. 



In Fig. [2j wc show that — at small e — the temperature 
scales develop as predicted from the multipole expan- 
sions. However, we note that the two-peak structure of 
the low-temperature specific heat persists all the way to 
the limit e — > 1 of short dipolar needles. This in par- 
ticular implies that the ground state of point like Ising 
dipoles on the kagome lattice with in-plane easy axes is 
given by the loop state. 

Magnetic monopoles in artificial spin ice: Given 
the efficient description provided by the monopoles for 
the phase diagram of the kagome array, it is natural 
to ask whether the same can be achieved for artificial 
(square) spin ice. Indeed, the first salient observation is 
that the ground state(s) of that array have zero charge 
on each vertex — this is possible as the local coordination 
is even, whereas it is odd for the kagome array. This 
makes contact with Wang's work, whose vertex types I 
and II correspond to Q = configurations. The ground 
state ("vacuum") thus contains no monopoles.^^ 

This leaves open the question whether monopoles 
nonetheless appear as effective low-energy degrees of free- 
dom, as they do in the three-dimensional spin ice mate- 
rials {Dy, Ho} 2Ti20T. As was pointed out in Ref. [l^, 
the short answer is no. Due to the 'antiferromagnetic' 
order in the ground state, monopoles cannot move inde- 
pendently as their separation leads to the creation of a 
costly domain wall: the monopoles remain confined. 

In our previous work^ however, we have proposed a 
modification of the square lattice geometry consisting of 
introducing a height offset between islands pointing along 
the two different lattice directions. If this offset, h, is 
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chosen so that the energies of all Q = vertices become 
degenerate — as e — >■ 0, the endpoints of the islands then 
form a tetrahedron, and h = — the ordering disap- 
pears, and the monopoles become free to move. Indeed, 
their interaction is given by 

V,ir) = -!^^^. (3) 
47r r 

Even though we have an array in d = 2. the interaction 
follows a c? = 3 Coulomb law, l/r, as the field lines are 
not confined to the plane. However, this is not yet quite 
the full story. In an appropriate low-temperature setting, 
there is in addition an entropic interaction — which hence 
scales with T — between the monopoles i^i^i This inter- 
action also is of Coulomb form, but it is two-dimensional 
in nature, as the entropic interactions due to the fiuc- 
tuations of the spin background do not know about the 
existence of a transverse third dimension: 

14(r) (xrin(^) . (4) 

The strength of such a confining potential can obviously 
be tuned by changing the temperature — at any rate, this 
potential is only weak and the monopoles will look essen- 
tially free on reasonable lengthscales. 

In experiment, this is true only in principle: a 
slowing-down of the dynamics makes the attainment of 
an equilibrium low-T state difiicultjii^ Indeed, the equili- 
bration dynamics of such arrays generally becomes very 
sluggish when one encounters energy barriers to fiipping 
individual islands. In the case of the kagome array, al- 
ready attaining the K2 charge-ordered crystal presents 



a formidable challenge in this sense. Knowledge of the 
states to be expected and their order parameters should 
nonetheless be useful in identifying equilibration strate- 
gies, especially since a magnetic field is a versatile probe 
of the spin ices due to their non-collinear easy axes<^ 

The scheme for describing dipolar arrays in terms of 
magnetic monopoles, as outlined above, obviously needs 
to be modified when ae becomes comparable to the trans- 
verse width of the needles, set to zero at the outset. How- 
ever, this will have a substantial effect only on the 0{l/e) 
term in the multipole expansion and it will in par- 
ticular not qualitatively change the relative sizes of the 
terms in the expansion in this limit. 

Finally, note that the systems studied experimen- 
tally typically have a linear extent of about a few hun- 
dred islands. Given the relatively long-range of the 
interaction — especially that a boundary need not be 
(magnetic) charge neutral — this means that even the 
ground state may contain domain structures interesting 
in their own right. 

Note added in proof. Recent work by Chern et ali^ on 
the kagome array obtains a phase diagram analogous to 
ours. 
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